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We clarify the questions rised by a recent example of a lattice Dirac operator 



l/-) , found by Chiu. We show that this operator belongs to a class based on the Cayley 

transformation and that this class on the finite lattice generally does not admit 
■ a nonvanishing index, while in the continuum limit, due to operator properties in 

Hilbert space, this defect is no longer there. Analogous observations are made for 
the chiral anomaly. We also elaborate on various aspects of the underlying sum rule 
for the index. 



1. Introduction 

Recently Chiu |IJ reported about a lattice Dirac operator, obeying the Ginsparg- Wilson 
(GW) relation @ and a mild condition on gauge field configurations, which does not admit 
a nonvanishing index, while everything else, in particular also the chiral anomaly, appears 
correct. The derivation of this was based on the sum rule || for the index. This remarkable 
observation rises generally the questions of what might be missing in the concept and how 
the apparent paradox resolves. To clarify these issues is the aim of the present paper. 

In Section 2, inspecting possible forms of unitary operators V which define GW Dirac 
operators, we note that the one implicit in Ref. [0J is a Cayley transform. In Section 3 we 
consider the general sum rule for the index and derive various relations which we need in 
the following. In Section 4 we establish the validity of the sum rule also in Hilbert space. 
In Section 5 we point out the differences to the Atiyah-Singer framework. In Sections 6 we 
generally analyze Cayley-type forms of V on the finite lattice as well as in the limit and 
completely clarify the issues in question. Section 7 contains conclusions and summary. 
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2. Forms of GW fermions 



For lattice Dirac operators D which satisfy the GW relation 

{^D} = p- 1 D To D (2.1) 

and are 75-hermitian 

£>t = 75J D 75 (2.2) 

due to O the remaining content of Q is p(D + Z>t) = DW = DDI The latter 
relation is the condition for unitarity of the operator V = 11 — p~ x D^ which is also 75- 
hermitian because D is. Therefore instead of using ( [2.1Q and ( |2.2| ) we equivalently require 
D to have the general form 

D = p(l — V) with V ] = V- 1 = 75 1/75 (2.3) 

and a real constant p. Since V is unitary, the spectrum of D is on a circle through zero 
around p. Real eigenvalues then are possible at and at 2 p. Clearly D is normal. 

There are three standard constructions of unitary operators V. They are up to constant 
phase factors given 

i) by X(\/X^X normalizing an operator X, not requiring particular mathematical 
properties of X (apart from X^X 7^ 0), 

ii) by the Cayley transform (Y — il)(Y + ill)" 1 of a hermitian operator Y, 

iii) by exp(iZ) with a hermitian generator Z. 

In addition to unitarity one gets 75-hermiticity of the constructed operators V by requiring 
X, iY, iZ to be 75-hermitian. 
Construction i) is used by Neuberger in the overlap Dirac operator |4[] choosing 



V = -X(VX^X)- 1 with X = ml + D w , -2 < m < , (2.4) 

where D\y is the massless Wilson-Dirac operator and where only X with X'X 7^ are 
admitted. This form is well known to give the correct results and does not need to be 
discussed further here. 

With respect to Construction ii) we observe that it is actually the basis of what has 
been introduced by Chiu |], [5|. In fact, the form D = a" 1 D c (l + rD c )~ 1 in Ref. |]], with 
an appropriate antihermitian operator D c and a suitable positive constant r , amounts to 
putting 

V = -{Y - it)(Y + it)- 1 with Y = irD c , (2.5) 

where p" 1 = 2ar relates the notations here and there. After a number of necessary 
preparations, we shall analyze Construction ii) in Section 6 in a general way and clarify 
the related questions. 
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Construction iii) so far has not been applied in the present context. For this one could 
use an operator Z of the form of Y mentioned above and tune r to cover the spectrum 
appropriately. There remain, however, problems with the periodicity to be settled. We 
will not address this case here. 



3. Index and spectrum 

To obtain the relations needed on the finite lattice, where D operates in a unitary space 
of finite dimension, we have the spectral representation D = J2k AfcPfc of normal D with 
eigenvalues A/, and orthogonal projections P^ = P\ . Using the 75-hermiticity of D, it gets 
the more detailed form 

D= E A J (if ) + /f ) )+ £ (A.Pf + A*Pj n) ) , (3.1) 

j (ImA i= 0) j (ImAj>0) 

satisfy = P,- 75 = ±P,^ and 75PP = P 7 - n) 75, one gets 



and one has P| + Pi projecting on the eigenspace with Ai = 0. Since the projections 
° (±) = P, (±) 75 = ±Pf } and 75P/ 



Tr(75P, (±) ) = ±N±(\ j ) , Tr( 75 P ; (l) ) = TV( 75 Pf ) ) = . (3.2) 

where the trace Tr is in full space and N + (Xj) and N_(Xj) denote the dimensions of the 
right-handed and of the left-handed eigenspace for real eigenvalue Xj, respectively. 

Using (|3.1|) with ( |3.2D the index A^ + (0) — A^_(0) is obtained from the resolvent (D — £ 1) _1 
by 

limTr(7 5 (-C)p-CIl)" 1 ) = iV+(0) - JV_ (0) (3.3) 

and one also finds 

limTr(7 5J D(P> - C l)" 1 ) = E - N -(^)) ■ (3-4) 

^ Aj^O real 

Adding up Eqs. ( |3 .3|) and (|3.4| ) the sum of their l.h.s. gets Tr(7 5 l) = and one obtains 



= N+(0) - M(0) + E (N+M - N -(*j)) ■ (3-5) 

Aj^O real 

This shows that to make a nonvanishing index at all possible, one has to require that the 
spectrum of D in addition to zero allows for at least one further real eigenvalue. 
The sum rule found by Chiu in Ref. B and used in Ref. M is the special case of (|3.5|), 



N+(0) - M(0) + N + {2p) - N_{2p) = , (3.6) 

which holds for Dirac operators satisfying the GW relation (|2.1|) . In that case (|3.4j) 
simplifies to 

(2p)- 1 Tr( 75J D) = N + (2p) - iV_(2p) . (3.7) 
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In the following it will be more appropriate to have the relations in terms of V, in which 
case p drops out completely. For this we note that by ( |2.3|) the eigenvalues and 2p of D 
correspond to the eigenvalues 1 and —1 of V, respectively. Thus, denoting the dimensions 
of the right-handed and the left-handed eigenspace at eigenvalue ±1 of V by iV + (±l) and 
iV_(±l), respectively, the sum rule ( |3.6| ) turns into 



N+(l) - iV_(l) + N+{-l) - iV-(-l) = , (3.8) 

requiring that in addition to +1 the eigenvalue — 1 of V can occur. Further, instead of 
fl3.7| ) with (5TB), we then have for the index 



N+{1) - JV_(1) = -N+(-l) + JVL(-l) = ^Tr( 75 ^) . (3.9) 
In this context it is instructive to note that the spectral representation of 75^ is 

75^ = Pi +) ~ Pt ] ~ Pt ] + Pt ] + E(H +] - Pt ] ) (3.10) 

k 

where and P 2 project onto the eigenspaces of V with eigenvalue +1 and — 1, 

respectively, and where p] ±] = ±(P fe (l) + P fc (n) ± e^ fc 7 5 P fc (l) ± e-^ fe 7 5 P fc (n) ) with < <p k < n 
(and [±] being related to signs of eigenvalues of 75 V). 



4. Sum rule in Hilbert space 

For the clarification of the questions raised in Ref. QX[ we have to make sure that the 
sum rule still holds in the continuum limit. Then the operators act in Hilbert space (as 
will be discussed in detail in Section 6), so that we firstly have to give a definition of 
the traces Tr in infinite space, too, and secondly to show that a continuous part of the 
spectrum does not contribute to the sum rule. 

The immediate prescription for the Tr-expressions is to consider them as the limit of the 
respective finite-lattice results. Technically more convenient is summing or integrating up 
a spectral representation after taking Tr of its contributions. Clearly any regularization 
can be used by which one arrives at the limit of the finite-lattice results. 

The Tr-expressions here typically involve differences of operators projecting on right- 
handed and left-handed space, respectively. This suggests to use the notion of cardinali- 
ties, which extends the one of dimensions to the infinite case. It allows to compare infinite 
spaces directly, using that two sets have the same cardinality if there is a bijective map- 
ping between them. For our purpose we restrict such mappings to ones between subspaces 
which have different chirality and call the resulting concept specialized cardinality. 

A simple example is Tr^l) = 0, which we have used to derive ( |3.5| ). In the infinite case 
one can decompose the identity operator as 11 = J2j Pj 5 with any orthogonal projections Pj 



4 



having finite trace Tr and being trivial in Dirac space, and define Tr^l) =: J2j Tr^Pj) = 
0. The special choice P n = 1 ®p(n), where p(n) is the projection related to lattice site 
n and lo the identity operator in Dirac and gauge field space, allows to make contact to 
the finite-lattice results. With respect to the general meaning, for 75 1 = P R — P L with 
the projections P R = |(1 + 75)! and P L = |(1 — 75)! a bijective mapping is provided 
by T = 74I, which gives T^P H T = P L and T^P L T = P R . This indicates that P R and P L 
project on subspaces which have the same specialized cardinality 

A further example is ^Tr^V) in ( |3.9|) . The discrete part of the spectral representation 
of the operator 75^ is given by (|3.10|) , while the continuous part does not contribute 
(as will be shown below). One notes that the subspaces on which J2kP^ and J2kPk~^ 
project have the same specialized cardinality. This follows because the projections p| + ' 
and pjp' are one-to-one associated and the subspaces on which such pairs project have 
both dimension |Tr(P^ + P^ ) • Thus only the part with P^ and P 2 contributes. 

We next address the fact that the spectrum of D now can also have a continuous part. 
The discrete part of the spectral representation is as before. The continuous part in terms 
of V is represented by a Stieltjes integral, 

V m = J\^dE contp = £e^dE contp - £e-*dE m . v , (4.1) 

where the subdivision into two integrals is possible because the projector function E conlfi is 
purely continuous. With V ] = 75^75 and P f on(p = E conip one gets ^ h E miLip = -P con _^7 5 . 
From this, using that E conlfi >E conlfi = E con(p E contp i = E contp holds for ip < <p', it follows that 
Tr(7 5 P con(p ) = . Therefore we obtain for any function P(V^. on ) 

Tr( 75 P(K con )) = f/(e^) d(Tr( 75 P con¥ ,)) = (4.2) 

and an analogous relation for the continuous part of D. Because the sum rule is based on 
expressions of type Tr^7 5 P(^)^ or Tr(^y 5 F(D)^ } it is indeed not affected by the occurrence 
of a continuous spectrum. 

5. Differences to the Atiyah-Singer case 

The sum rule, which in Ref. |j] and in the present paper is seen to be of crucial 
importance, exhibits the mechanism allowing for a nonvanishing index on the lattice. 
This mechanism is related to the chiral noninvariance, which is conceptually necessary 
in quantum field theory, and in its lattice definition is present from the start. Clearly in 
lattice theory the space structure itself neither does depend on the Dirac operator nor 
does it get chirally asymmetric. 

This is in contrast to what one has in case of the Atiyah-Singer Dirac operator , where 
to admit a nonvanishing index the space structure itself must be allowed to depend on the 



5 



particular gauge field configuration. Because this appears to be not sufficiently realized, 
we briefly point out some details. 

The definition in the Atiyah-Singer case is based on (Weyl) operators P A s and D A ~s 
which map from the total right-handed space 8^ to the total left-handed space £^ and 
back, respectively. It is given in the combined space by P AS = D^+bQ where 

D$ = on £W and D A f = on £&\ Because of P A s )f = the Dirac operator 
P AS is self-adjoint and since it acts on a compact manifold its spectrum is discrete. Thus 
it is represented by P AS = J2j A AS jP AS j which implies D 2 AS = J2j ^isj^sj- 

On the other hand, one gets D AS = Das*D a ~^ + -Pas Pas , which can be evaluated not- 
ing that the operators P AS ' ) P A s' ) and P A t' ) P A s' ) map within £^ and £^~\ respectively, 
and are selfadjoint and nonnegative. With the eigenequation P(~)p( + )<I>j = Kj<&j in 

£(+) 

one gets the eigenequation pWp( )(P( + )$j) = Kj(D^ + ^j) in £( '. Further, from 
(p (+) $j r / 1 p(+)$ jr ) = ($ ir /|_D ( ~ ) L> (+) $ J v) = Kj($j r i\$j r ) one sees that, except for Kj = 0, 
for each of the common eigenvalues Kj the eigenspaces must have the same dimension. 
Therefore, except for Kj = 0, the operators Pas Pas and Pas Pas have the same spec- 
tra. Denoting the projections on their eigenspaces by P^j and P| s j, respectively, and 
comparing the above expressions for D AS it then follows that D AS = J2j ^as jC^is] + ^Is]) 
and that one always has 

iV + (A ASi ) = i\L(A AS i) for A ASi ^0 (5.1) 

for the dimensions N + (\ AS j) = TrPf^j and AL(A AS j ) = TrP| s j of the eigenspaces. 

For Kj = \ Asj = one gets kerP A t } = kerPis ) P A t ) and kerP A s } = kerP A t ) Pis ) 
(as is obvious from left to right and follows from right to left from ($|P A g' ) P A t' ) $) = 
{DaV^IDaV®)). Thus the eigenspaces with eigenvalue zero have the dimensions iV±(0) = 
dim kerP A ^ and the index becomes N+(0) — N-(0) = dim ker D A V — dim ker Pas^- 

It is now seen that because of (|5.1|) , to admit a nonvanishing index iV + (0) — iV_(0), the 
space structure itself must be allowed to vary and to get chirally asymmetric, and thus to 
depend on the particular gauge field configuration. It also becomes obvious that in case of 
a nonvanishing index the specialized cardinalities of £^ and £^ get different. Indeed, 
while right-handed and left-handed eigenspaces with the same dimension are associated 
at each A AS j ^ 0, this then does not hold for those with A ASJ - = 0. 

In addition to the fundamentally different space structure, a further difference is that in 
the Atiyah-Singer framework one considers a differential operator on a compact manifold, 
while in the definition of quantized theory one deals with a lattice operator and its subtle 
limit on R 4 . Thus altogether special care appears appropriate when considering analogies. 

6. Analysis of Cay ley-type V 

To investigate Construction ii) of Section 2 we consider the general choice 

V = -{Y - it)(Y + it)- 1 = 2(1 + Y 2 Y l - t + i2Y(t + Y 2 y l . (6.1) 
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In (|6.1|) for the eigenvalue y = of Y one obviously gets the eigenvalue v = 1 of V. 

On the finite lattice the operators act in a unitary space of finite dimension. Therefore, 
requiring Y to be a well-defined hermitian operator, its spectrum consists of a finite num- 
ber of real eigenvalues (which are discrete and finite). Introducing s = max(|y min |, |y max |), 
where y min and y max denote the smallest and the largest eigenvalue of Y, respectively, ( J6.1|) 
gives 

2 2s 

Re v > - — 1 for all v , lim v\ > for Re v < . (6-2) 

1 + s 2 1 + s 2 

We see from this that the eigenvalue v = —1 of V cannot be reached. Thus on the finite 
lattice Construction ii) does generally not meet the basic requirement needed according 
to ( p.8|) to allow for a nonvanishing index. 

In addition, we observe that , because of not reaching v = —1, one always has P 2 = 
in ( |3.10| ). Thus on the finite lattice already the chiral anomaly turns out to be affected. 



The obvious obstacle which prevents from reaching the eigenvalue — 1 of V is that on the 
finite lattice Y is bounded. A related problem is that there the inverse Cayley transform, 

Y = —i(V — t)(V + I)' 1 , (6.3) 

is not valid for all unitary operators but only for the subset for which the spectrum does 
not extend to — 1 . 

The crucial observation now is that the indicated restrictions no longer hold in Hilbert 
space, where one gets a well-defined connection between general unitary operators V and 
selfadjoint operators Y which can also be unbounded. To recall how this comes about we 
start from the general spectral representation of unitary operators, 



V = / e i,p dE v , (6.4) 

where the projection function accounts for discrete as well as for continuous con- 
tributions. Naive insertion of ( |6.4|) into ( |6.3|) does not generally make sense because 
— i(e tv — l)(e lLp + 1) _1 = tan| is not bounded, diverging for = ±tt where the value — 1 
of the spectrum of V is reached. However, Y is well-defined on Hilbert space vectors / by 

Yf= lim f tan^d(iV) (6-5) 

in the sense of strong convergence. This is seen noting that with / = (1 + V)g one gets 
tan 2 ^ d| \E<pf\ | 2 = 4 sin 2 ^ d| l-E^I | over which the integral from — it to n is obviously finite. 



Thus with unbounded operators Y in (6.5) we indeed get unitary operators V in (|6.4j) 



with a spectrum extending to — 1 , as is necessary for proper working of the sum rule 
( |3.8|) and to guarantee correct results in ( |3.9|) . It is seen that for this a Hilbert space 
is necessary, which not only has infinite dimension but also includes its limit elements. 
Obviously the latter here is of crucial importance. 
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Since the projection function E v is common to (|6.4j) and ( |6.5D , we can trace back for V 
as well as for Y to the projections and P 2 (related to the eigenvalues +1 and —1 of 
V, respectively) which only enter the sum rule. With E v being strongly continuous from 
the right, lim e ^ E v+0 = E v , for this one has the relations \im E ^ (E — E_ £ ) = P^ +P[ 
and lim e ^ (^ - E^ e ) = P 2 {+) + P 2 H '. 

To see how the Hilbert space of interest arises in the continuum limit, one has to 
reconsider things from the appropriate point of view : On the infinite lattice one gets 
the Hilbert space of sequences 1 2 (in a basis of which a vector is related to a lattice 
site, a Dirac index and a gauge field index). The unitarily equivalent Hilbert space 
L 2 (7r; k) of functions f{k) with — n < < n (and Dirac and gauge-group indices being 
suppressed) is obtained from I2 by a Fourier transformation. Introducing the lattice 
spacing a and variables p = k/a the space L 2 (7r;fc) becomes L 2 (7r/a;p). By the limit 
a — ► one then gets the operators in L 2 (oo;p) from the ones in L 2 (7r/a;p). The space 
L 2 (oo; x), unitarily equivalent to L 2 (oo;p), is again obtained by a Fourier transformation. 
Instead of proceeding in the more instructive way sketched, one can also realize the direct 
way from 1 2 to L 2 (oo; x). The equivalent spaces L 2 (oo; x) and L 2 (oo;p) are the ones one 
has in the continuum limit. 

In detail the definition of the operators of interest by the indicated limit needs some 
care. Firstly, since the limit element cannot be given explicitly, we resort to the definition 
by all matrix elements, i.e. by weak operator convergence. Secondly, because two spaces 
are involved, this weak limit is to be slightly generalized. To show that this can be 
properly done, we introduce f a (p) = /(p)II M 0(7r/a — p^) with f(p) e L 2 (oo;p) and the 
operator O a , requiring O a (p',p) of L 2 (oo;p) to be equal to O a (p',p) of L 2 (7r/a;p) for 
-a /a < Pf, < 7i /a. Then (f a \O a g a ) in L 2 (oo;p) equals (f a \O a g a ) in L 2 (7i/a;p) for all 
finite a and (f a \O a g a ) — > (f\Og) for a — > defines the operator O in L 2 (oo;p). 

Thus, though we cannot give explicit expressions for the operators in L 2 (oo;p) and 
L 2 (oo;x), we can unambiguously deal with them: First calculating the desired matrix 
elements of the operators or of the functions of operators of interest on the infinite lattice 
and then performing the a — > limit leads to the results in terms of the respective matrix 
elements. It is seen here that what is practice in correct approaches can be precisely 
formulated and justified in Hilbert space. 

The mapping from O a to O is not invertible and the spectra of O a and O can be 
substantially different. Obviously the momenta get unbounded in the limit, and this also 
holds for the operators Y of physical interest. Therefore by the connection between ( |6.5| ) 
and flO]) the spectrum of V can extend to — 1 , as has been shown to be necessary for 
correct results. 

We now see how the paradox of Ref. Jl| resolves : The chiral anomaly there, being 
already a continuum quantity, should not be considered with the sum rule on the finite 
lattice but with the one valid in the continuum limit. In this limit the chiral anomaly as 
well as the sum rule can get correct properties, while on the finite lattice both of them 
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are affected by not reaching the eigenvalue —1 of V . The problems on the finite lattice 
have been seen to occur not only in the example of Ref. |J but for any well-defined Y. 

7. Conclusions and summary 

Inspecting possible forms of unitary operators V which define GW Dirac operators, we 
have pointed out that the Dirac operator proposed by Chiu is based on an operator V of 
the Cayley-transform type. After a number of necessary preparations, we have analyzed 
the properties of this class of operators. 

Considering the general sum rule for the index we have derived various relations in finite 
space which have been important in the following. To make sure that the sum rule still 
holds in the continuum limit, we firstly have shown that the extension to infinite space 
does not spoil it and secondly that a continuous part of the spectrum does not contribute 
to it. 

We have stressed that the mechanism for admitting a nonvanishing index, exhibited by 
the sum rule here, is different from that in the Atiyah-Singer case, where instead the space 
structure itself is allowed to depend on the particular gauge field configuration. 

Considering Cayley-type operators V we have seen that for them on the finite lattice the 
sum rule generally does not allow for a nonvanishing index (and not only in the example 
of Chiu) because the eigenvalue — 1 of V cannot be reached. In addition, we have seen 
that in case of such V on the finite lattice already the chiral anomaly is affected by this. 

Turning to Hilbert space we have pointed out that there the Cayley transformation 
provides a well-defined connection between unbounded self-adjoint operators and general 
unitary operators V for which the eigenvalue —1 can be reached. We then have shown 
that, after introducing a slightly generalized weak operator convergence, the continuum 
limit precisely provides the appropriate Hilbert-space formulation. 

Our results make clear that in the continuum limit the problems of the sum rule for the 
index and of the chiral anomaly with Cayley-type operators V disappear. The paradox 
of Ref. |l[ has been seen to resolve in that the chiral anomaly there, being already a 
continuum quantity, should not be considered with the sum rule on the finite lattice but 
with the one valid in the limit. 
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